arXiv: 1501.0774lv2 [math-ph] 12 Mar 2015 


NEW PERIODIC SOLUTIONS FOR NEWTONIAN n-BODY 
PROBLEMS WITH DIHEDRAL GROUP SYMMETRY AND 
TOPOLOGICAL CONSTRAINTS 


ZHIQIANG WANG AND SHIQING ZHANG 

Abstract. In this paper, we prove the existence of a family of new non¬ 
collision periodic solutions for the classical Newtonian n-body problems. 
In our assumption, the n = 21 > 4 particles are invariant under the di¬ 
hedral rotation group D; in R 3 such that, at each instant, the n particles 
form two twisted 1-regular polygons. Our approach is variational min¬ 
imizing method and we show that the minimizers are collision-free by 
level estimates and local deformations. 

Keywords: Periodic Solutions, Newtonian n-body Problems, 
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1. Introduction and Main Result 

Many authors (for example [6] [8] [9] [18] [22j ) used the variational method 
to discover many new periodic solutions of the classical Newtonian n-body 
problems in the last fifteen years. In particular, Chenciner and Montgomery 
[6] proved the existence of the remarkable figure-8 type periodic solution for 
planar Newtonian 3-body problems with equal masses. Ferrario and Ter- 
racini [8] simplified and developed Marchal’s m important works and in¬ 
troduced the rotating circle property, proved that if the motion has certain 
symmetry under some group action having the rotation circle property, the 
solution exists and has no collision. Also Fusco et al.|9] proved the exis¬ 
tence and collisionless of a number of new and interesting motions with the 
invariance of certain platonic polyhedra group action and some topological 
constraints. 

In this paper, we consider a system of n = 21 positive masses with 
their positions x(t ) = (xi(t),X 2 (t),... ,x n (t)) T moving in the space under 
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Newton’s law of gravitation: 


( 1 . 1 ) 


M.x(t ) = 


dU(x(t)) 

dx 


where the potential function U(x) = Y , and Ad = diag{mi,rri 2 , ■ ■ ■, }. 

It is well known that looking for periodic solutions for (11.11) is equivalent 
to seeking the critical points of the Lagrange functional A : A —>• RU {+oo} 


A(x(t)) = / C(x)dt = / (K + U)dt 


rT n 


E 


1 


rmiij 


'0 


i =1 


+ E 

Kj 


miirij 
\xj — Xj 


■)dt 


on the set 

A = { x {t) € ^(R/TZ.R 3 ")!®^*) ± Xj(t),Vi ± j,Vt € 
and our approach is based on the following basic lemma: 


Lemma 1.1. f|17j) Let X be a reflexive Banach space, M C X is a weakly 
closed subset, f : M —>• R is weakly lower semi-continuous; if f is coercive, 
that is, f(x) —>• +oo as |x| —> +oo, then f attains its infimum on M. 


There are two difficulties in this approach: one is the lack of coercivity 
of the action functional A on the whole set A; and the other is that in critical 
points, there might be trajectories with collisions. To obtain the coercivity, 
one can consider the functional A on some symmetric subspace Ag> C A 
such that _4 .|a g is coercive. And the following famous lemma proves that 
the critical point on Aq is also a critical point on the whole space A. 


Lemma 1.2. (Palais principle of symmetric criticality pAj ) 

Let G be an orthogonal group on a Hilbert space A. Define the fixed 
point space: Aq = {x € A\g ■ x = x, \/g E G }; if f E C ,1 (A, R ) and satisfies 
f(g-x) = f(x) for any g € G and x € A, then the critical point of f restricted 
on Ag is also a critical point of f on A. 


Here is a traditional way to define the group action on loop space A 
such that A(g-x) = A(x). Let G be a finite group with three representations 
p : G —» 0(d), t : G —» 0(2) and a : G —>• such that 

g-x{t) = {p{g)x (j(g -i m {T{g- 1 )t),..., p{g)x a{g -i )[n) {T{g- 1 )t)). 
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Where we only consider homomorphisms a with property that Vg E G : 
(c 7 (g)( 1 ) = j =>- 7?ij = rrij) and for more detail, we refer the readers to [ 8 J. 
Also we can add some topological constraints on A.q to get an open cone 
1C C A g with the property 


dIC C A g = {u E A g : 3t c € R, i / j : Ui(t c ) = r^(t c )}, 


and find critical points inside 1C ora usd. By the above arguments, we 
can distinguish geometrically different solutions, and get the coercivity even 
if vA|a g is not coercive. Indeed, if we are able to prove that a minimizer u * of 
-4| Tc exists and for any collision trajectories u c E dIC C Aq: A(u*) < A(u c ), 
then we must have it* E Ag and is collision free. Thus by Lemma [L2l u * E 1C 
is a critical point of *4 |a and therefore a solution of the n-body problem. 

Suppose the motions are in the space O — and let e,j be the unit 

vectors of the coordinate axes for j = 1,2,3. Denote the rotation of 

^cos — sin ^ 0 \ 


angle around /( 3 -axis R = 


around £i-axis S = 


sin cos 0 


0 


0 


, the rotation of angle 7 r 


1 / 



. Then Di = (R , S) is the dihedral group 


OH]) of order n = 21 and it is a group of rotations with their rotation axes 

klT 


r = £3 where L k is the line 


£2 = 6 tan—, 


>3=0 


£3 = 0 . 


For simplicity, we denote Di = {Rj}™ =0 , where R k = R k , Ri+k = R k S 
for k = 0 , 1 ,..., l — 1 , i.e. 



/cos^ 

-sin^ 

(A 


/cos^ 

sin^ 

° ^ 

Rk = 

sin 

cos^ 

0 

! + — 

sin 

-cos^ 

0 


V 0 

0 

1/ 


V 0 

0 

-1/ 


Now we consider n = 21 > 4 point particles u(t) = (uo(t), ui(t ),... , u n - 1 W) 
with equal masses in space with the following symmetry: 


(1.2) Uj(t) = RjU 0 (t), j = 0, l,...,n- 1. 

Under this symmetry the trajectories have the property that, at each instant, 
the n point particles form a two nested regular /-polygons with the same size. 
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By letting the the masses rrii = 1 and under the assumption of the symmetry 
m, the action functional can be written as 

A(u(t)) = A(Mt)) = (m*)I 2 + E 

In [7], Ferrario and Portaluri studied central configurations with this dihe¬ 
dral symmetry. In Section 2 of [9], Fusco et al. got a new periodic solution 
in the case of n = 4 by applying some topological constraint, here our result 
is a generalization of theirs. Obviously, the trajectories u(t) are uniquely 
determined by the trajectory uo(t), and in [9], uo(t) is called the generating 
particle of the motion. Also we need some other symmetric constraints on 
the loop of uq: 


(1.3) 


u 0 (t) = Rouo(-t), 

T 

u 0 (t) = R s u 0 (— - t). 


Where s,h < l are some positive integers and R & = sin — cos 


kn 
l > 


/cos^ 

sin^p 

sin^ 

-cos^ 

V 0 

0 

A; = 0,1,... 

,1-1. 

>(* - D = 

R s uo{t - 


T\ 

h> 


Since uo(t) is a T-periodic solution i.e. uo(t) = uo(t + T), we must have 
Rg = Rq = id which implies 


(1.4) 


sh = l mod l. 


Moreover if T is the minimal positive period, then h is the minimal positive 
integer satisfying (|1.4I) and of course h\l. Also applying (|1.2I) we see that 
u 0 (t) = u s (t - f) and 

ui(t) = Suo(t) = SR s u 0 (t - ^-) = Ri_ s Su 0 (t - ^) 

T T 

= R 2 l- s Uo{t — —) = u 2 i-s(t - —). 

h h 
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If we denote j to be the least nonnegative residue of j modulo l , i.e. j = j 
mod l and 0 < j < l, we have 

Mt) = u s(t - T ) = U Ts {t -—) = ■■■= 

Mt) = ^i+T(* - |) = %i+r(t - = " ■ = u {h _ 1)s+1 (t - ^ h ^ T ), 

T 2 T (h — 1 )T 

u i/h- 1(0 = — = ^s+Z/Ti-l^ ~h^ = "* = u (h-i)s+l/h -^ )’ 

T 2T (h— 1)T 

< ui(t) = u 2 i-s(t --j;)= u 2 i-Tstt - ~fr) = • • • = - 1 -)> 

T 2 T (h — 1)T 

u 2l-l(t) = “2Z-(i+I)(* “ = U 2/-(2i+l)(^ ^") = ■ ■ ■ = U 2«-((/i-1)8+1)^ h )> 


u 2l-l/h+l(t) — U 2Z-(s+Z//i-l)^ ~ Un 


2l-(2s+l/h-l) 

tx {h-\)T s 

~ U 2l-((h-l)s+l/h-l)V ; )' 


, 2 T 


That is to say the n = 21 particles’ motion are composed of 21/h choreog¬ 
raphy trajectories. For example when l = 12, s = 9, then we have h = 4 
and 


T T 3T 

«o(t) = «9(* - j) = Mt - = u 3 (t - —), 

T T 3 T 

Mt ) = «io(* - j) = «7(* - = Mt - 

T T 3 T 

Mt) = u n(t - -) = u 8 (t - ^0 = «5(* - 

T T 3T 

u 12 {t) = ui 5 (t - —) = ni 8 (t - —) = « 2 i(i - —), 

r T 3 T 

U23 (t) = U 14 (t - —) = Rl7(i - 2-) = ^2o(t - —), 

T T 3T 

U22 (t) = Ui 3 (t - —) = Ui 6 (t - —) = Ui 9 (f - —). 


Remark 1.3. We notice that the angle between the plane P s and Po is the 
same with the angle between Pi- S and Po> and hs = l mod l if and only if 
(l — s)h = l mod l which means l — s and s implies the same h. So in (|1.3I) . 
the two symmetric constraints are the same in geometry, and we can only 
consider the integer 1 < s < ^ = n/4. 
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Remark 1.4. By the condition (11.31) . we see that I = [0, is a fundamen¬ 
tal domain of dihedral type for the trajectories (see 0 ) for more details), 
which implies that the motion of the particles on the whole period [0, T] is 
determined by their motion on I = [ 0 , through the symmetric conditions. 
And from ( 11 . 21 ) we see that uq is the generating particle, so in Section [3] we 
can only consider the motion of the generating particle Uq in the interval 

I = ^11- 

Let G s = (R, S, R s , Ro) with the following representations: 
p(R) = R , r(R)t = t, a(R) = (0,1,..., l — 1)(Z, l + 1,..., n — 1), 

i -1 

p{S) = S, r(S)t = t, a(S) = (0, l) n (I - M + k), 

k =1 

p(R s ) = R s ,r(R s )t = ^ - t, cr(R s ) = id, 
h 

p(Ro) = Ro, r{Ro)t = -t, a(R 0 ) = id, 
and set 

A s = {u(t) € A : u(t) satisfy ( 11 . 21 ) (II. 3p }. 

It is easy to check that A s = A q s , i.e. looking for trajectories with properties 
( 11 . 21 ) (11.31) is equivalent to seeking for critical point of A on Ag s . 

In 0, Ferrario and Terracini proved that *4 .|a g is coercive if and only 
if Xq = {x € X\g ■ x = x,Mg £ G} = 0 where X is the configuration space 
of the particles. Obviously, in our assumption, ^4 |a 0 is not coercive since 
(e 3 ,..., e 3 , —e 3 ,..., — es) € Xq where e 3 = (0,0,1). So motivated by [9j, we 
add some topological condition on A g s , to get the open cone IC S described 
in the previous. From (11.31) we see that rto(0) € Po and uo(^) € P s , so we 
let 

(1.5) JC S = { u(t) € A s : u 0 ( 0) G P 0 ".' u o(^) € P+} 

where we have set P 0 _ = {p € Po : p ■ e% < 0} and P+ = {p € P s : p ■ e% > 0}. 
Now we state our main theorem: 

Theorem 1.5. Forn = 21 > 4 and every integer s < max{l, () 3//2 ^72 log 
1} where 7 « 0.57721566490153286 is the Euler-Mascheroni constant, there 
exists a T-periodic solution u * € IC S of the classical n-body problem. 
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Corollary 1.6. For n = 21 > 4 and s = 1, then h = l and there exists a 
T-periodic solution u * E K,\ of the classical Newtonian n-body problem. 

The case for n = 4 was discussed in Section 2 of |9], here we generalize 
their result. By Remark 1 1.31 s can be chosen larger than 1 when n > 8. If we 
let f(n) = 3/2 ^72 log n +7 “ we see tliat /( n ) is monotone increasing 

and we have /(4) » 0.4697, /( 6 ) « 1.1400, /( 8 ) ss 1.7376, /(10) ss 2.2931, 
/(14) « 3.3262, /(26) ~ 6.0995, etc. That is to say, for n > 10, we can 
choose s = 2 such that there exists a T-periodic solution w* E 1C 2 of the 
classical n-body problem. 

Remark 1.7. Actually, the value of integer s depends on the estimate of 
excluding total collisions. By doing more explicit computation we can prove 
that, for n = 8, there exists a T-periodic solution 74 * E /C 2 of the classical 
8-body problem. 


2. Coercivity 


Proposition 2.1. A\jc s is coercive and dlC s C A q s . 

Proof. From (11.21) we see that there is a collision if and only if there exists 
some t c E R such that Uo(t c ) E T, where T are the rotating axes of the 
dihedral group Di. So it is obvious that d)C s C A g s and in the following we 
will prove the coercivity. 

Applied Newton-Leibniz Formula and Holder Inequality, we have 



T 


T 


T 


< (T)kjf“ Mill 2 *)* 


Since uq(0) E P 0 ,uq(^)€ we must have 
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which induces that 


\\Mt)\\h > C^-) - u o (0)\ 2 

T 


>Cl 

— C 1 l 

> C\u 0 (0)\ 2 , 


T 


SIT 


M^)| + M°)I - 2|* io (0)|M—)|cos 


oh r T "1 2 QHT 

I Wo (0)1 COS — - |it 0 (—)| + Ci|w 0 (0)| 2 sin 2 — 


which implies that H'i^oHz, 2 is an equivalent norm of H 1 = W 1,2 and the 
coercivity for the functional A follows. □ 


3. Estimate on Collisions and the Proof of Theorem 11.51 

In this section, we show that the minimizer it* G K s is free of collisions. 
Firstly we exclude total collisions and in Section f3.2l we discuss the partial 
collisions. 

3.1. Total Collision. Our way to show that there is no total collision is 
based on level estimate. Firstly assuming that a total collision happens 
in ti* G JC S , we show that there is a lower bound of the action functional 
B < A(u *), and we construct a test loop u € K s without collisions such that 
A(u) < B. Then we have A(u) < A(u*) which contradicts with that it* is a 
minimizer. Thus the minimizers u* G K s is free of total collisions. 

Lemma 3.1. (Gordon’s Theorem pTjJ Letx G W 1,2 ([ii, £ 2 ], R d ) andx{t\) = 
x{t 2 ) = 0. Then for any a > 0, we have 

f t2 ,1, . ,n a 3 . . 2 2 1 

/ (ol^l 2 + rr)^ — o( 27r ) 3a3 ( i 2 - ti)s. 

J t\ 2 \x\ 2 

Proposition 3.2. (Lower bound estimates for A with total collisions) 
Assume that u G K s has a total collision. Then 

^(»)> 3 ( ’‘ 4 ~ 1 ) (2ft)LLIrU B . 

Proof. Since the center of mass is at origin, the functional A can be written 
as 
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where the total mass M = ^mj. This formulation came from ([19]|20|) 
and has been widely used to obtain the lower bound estimate of collision 
paths ([3].[5].[22].|21]etch In our assumption, we have m\ = • • • = m n = 1 
and M = n. Moreover, if u £ IC S has a total collision, then they collide at 
least h times in the interval [0, T). Applying Lemma 13. 11 we have 


A(u) = 


1 


n 


u 


r l, 


n 


> 


i<j 

1 n(n — 1) 3 


\Uj — u- 


r\dt 


n 2 
3 (n - 1) 


2 2 T 1 

2 (2tt) ln§(-)5 x h 


(2/i)3vrin3T3 = B. 


□ 


Next we construct test loops u € K s such that A(u) < B in two different 
ways. The idea of the first one is from Fusco et al.[9], but it holds only for 
s = 1. The second one holds for s < () 3 ^ 2 ^72 i og n +7 ~ anc ^ ^ nee ds 
some more explicit analysis on the potential U. 

Proposition 3.3. (Upper bound estimate for s = 1) 

When s = 1, there exists u € JCi such that 

3 1 221 

A(u) < -(2 h 2 )^n(n — l) 37 r 3 T 3 . 

Proof. We construct the test loop u similar to that in Proposition 5.3 in 
[9j. Assume the generating particle uq moves with constant speed on a 
curve which is the union of two quarters of circumferences C 1 , C 2 of radius 
r tan ^. C\ has the center (r, 0,0) and lies on the plane = r. C 2 has 
the center (r cos f, r sin f, 0) and lies on the plane £icos f + £2 sin j = r, see 
figure 1. So the constant speed of generating particle uq = , and the 

kinetic energy 

*<4> = 

From the definition of u, we see that \u.i — Uj\ > 2r tan ^ for all i A j\ 
which implies 

n(n — 1) 


V(u) < 


1 

2rtan : 


therefore we have 


. . f T , . 2h 2 mr 2 r 2 tan 2 - n(n — 1 )T 

A(u) = / K(u) + V(u)dt < ---^ 1 j 

Jo 1 


4r tan - 
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Figure 1 . 


and the conclusion follows if we choose r 


( 16 h 2 )3 tan ^ir 2 / 3 


(n— 1 ) 1/,3 T 2 / 3 

—TT-~ 


□ 


From Proposition 13.21 and 13.31 we see that < { 2 (n-i) ) 1 ^ 3 < i- e - 
A(u) < B. Next we consider the situation s > 2, since s < j, in the 
following we suppose n = 21 > 8. 

Proposition 3.4. (Upper bound estimate for the general s) 

When n = 21 > 8 , there is a test loop u € IC S such that A{u) < B for 


s + 1<(^) 3/2 


2 3 / 2 log n+7 ' 


Proof. Assume that the particle uo moves with constant speed on the sphere 
|Sol = a, with the radius a to be determined later. More precisely, suppose 
uo(p(t),8(t)) = a(cosipe d ^^~^,siiup(t)) where 




(<p(t),0(t)) = (-- ,ut ) 
n 



(.<p(t),0(t)) = (ut 


2n 7T 


^ ^4 h(s + 1) ’ 4 h(s + 1) 
3 T T 

t £ Ah(s + 1)’ 2h 


T 3 T 


n n 


7T Z7T 


n 


n 


where u = = 2/i ^y 1 - >7r , see Figure 2. 

First we claim that, for every t E [0, X], the potential 


77, 

U(u(t)) < -—(logn + 7). 
zan 


(3.2) 
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£3 



And the kinetic energy 


im < 


2h 2 na 2 (s + 1) 2 7 t 2 

PT 2 


so the functional 

f T r , [ T 2h?na 2 (s + l) 2 n 2 n 2 

A ^ = Jo [ R+U } dt< J 0 [- pf 2-+ — (logn + 7)]dt. 

we choose a = VlUL-L (i°g»+7) 1 (l) 2 / 3 then 

2tt(s+ 1)3 n 

A{u) < ^(y)inl(logn + 7) 2/3 (s + l) 2 / 3 T 1/3 , 


A(u) 

B 


n 2 
n — 17T 2 / 3 


( 


log n + 7 
n 


) 2/3 (s + 1) 2/3 < 1, 


and the conclusion follows. 


□ 


Now we prove the estimate (13.21) . first we need some more explicit es¬ 
timate on the potential U. We notice that, at every instant, the n par¬ 
ticles form a twisted regular /-polygons by the symmetric conditions. Let 
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uq = (a cos a sin +), then the potential function can be written as 

(3.3) 


n— 1 


U = U (uq) = U(tp, 0) = ^ \( R j - i?o)tt 0 r 1 


1 =1 


l-l 


= 2Z( cos< ^ I 1 + X][ cos2< ^l 1 “ e 2(?xAT ^| 2 + 4 sin 2 (^] 1/2 } 

i=i i=i 

= i ^{20+g[ S m 2( lf -0) + tan+]"i}, 

where C/ = |1 - £/I _1 = h Y ! 3 ~=i csc T and & = 

Also we notice that 

cos 2 +|1 — e - 2 0\/-i£i| 2 _)_ 4 s i n 2 ^ 

= cos 2 +(2 - e- 26 ^ 1 ^ - e 26 ^ 1 ^) + 4 sin 2 p 

= 2 + 2 sin 2 <p - cos 2 ip(e~ 2e ' / ^ 1 ^ + e 26 ^^’) 

- (1 + m _ A_£ll+ (e -^ + e ^p )} 

(l + sm<^) 2 (l + sm<y9) 2 

= (1 + sin ^ [l + (+£++ - lz^ {e - 20 V=i ( i + e »l=^+ 

= (1 +sin v ) 2 |l - l+++- 2 »' /=T {+. 

1 + sm (p 


That is 


u(tp,e) = 


n 


2 a cos (p 


Q + 


COS (f 


En 


l-sin^ M/ rr 


1 + sin p +-;' 1 + sin 99 


ii ’ 


let r = ^ and £ = e 26, v / ^’ ; then 


(3.4) 


C/(r,0 = d^(C, + lF^|l-r+r'), 


da-y/r 


1=1 


since there is an integral representation of )Th =| |1 — r^j | 1 , we can state 

Lemma 3.5. For r € (0,1) and £ = e -20v/ ^" the potential U(r , £) can 6e 
written as 


(3.5) tf(r,£) = 


n 


( 1 + r ) /\/f f 1 (1 1/,2 £ ^ 2 1 — ( ir ) 


2Z 


Q + 


Aay/r V 7i J o (1 — tr 2 ) 1 / 2 |1 — (fr^yp 


dt 
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Proof. In 0, Ferrario and Portaluri gave a general integral representation 
of y'Sj—j |1 — with 0 < a < 2. Since this lemma is rather important 

to our analysis, we prove it here again in the case of a = 1. 


|l-r££/| 1 = (1 -rfff) 1/2 (1 - r£ 1 ? / J ) 1/2 



(-rfV) fc 


2 1 r k+h \ (ccP n 


mr 


E **(«?)"■ 


since T(^) = y/n, and the Beta function B(x, y) = t x 1 (l — t) y 1 dt = 

'(v) 

r (x+y) 


r(x)r (y) we } iave 


and 



k\ 

* r ( fc + ^) 


= (-i) 


r(fc + i)r(i) 
(-i)"r(Hi)r(i) 


bn = (-l) n E 


7r 

r(/c +1) 


r tk , 

7T 

/0 (l-t)W/2 

(-D* 

/•' 


dt, 


dt 


k—h=n 

k,h>0 


7T J o (1 - t)l/2 t l/2 l 


2 i 


£^1 

/i=0 U 


h /*1 ^/i+n 


(1 - i)l/2*l/2 l h 


2 1 .2h+n^ 


{tr) 7 


it J 0 (l-ty/H 1 / 2 


(1 - tr 2 )~ 1/2 dt, 
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since b n 


b-n and Y!j=i PP' = 


if n ^ 0 mod l, 

so we have 

if n = 0 mod l, 


1 

7 


i 


Eh 


r «fi = 7E E m{{/)” = E«'”+E« 


l oo 


—In 


n= 1 


j=l n=—oo n=0 


\ln 

, (1 - tr 2 )~ 1,2 C ln dt 


(tr 


n =0 

oo j 

+ S’ ? - / » (i-*) i/2 « i/2 ’ 

' 2 n=0 n=l 

1 r 1 (1 - i)- 1 /^- 1 / 2 1 - (tr) 2/ 

7T J 0 (1 — tr 2 ) 1 / 2 |1 — (£?’£)* | 2 


rdf. 


Thus the integral representation (13.511 holds. 

□ 


Lemma 3.6. The potential U of the test loop u 


U{u) < max{f/(0. J)> *7(^.0)}. 


Proof. By ()3.5I) 


U{r,6) = 


n( 


(1 + r) r ly/r f 1 (1 — t) l ^ 2 t 1 / 2 


Ci+- 


L 


1 — (tr) 


21 


\ayjr V 7r J 0 (1 — tr 2 ) 1 / 2 1 + (tr) 21 — 2 (tr) 1 cos(2 W) 


dt 


which implies ^ < 0 for 0 < 0 < Jfj, 0 < p < Thus we get for given 

0 < p < §, 


sup U(p, 6) = U(<p, 0). 

0e[o,§l 


From (13.311 . we have 


3C7 

dip 


n sm ip 
2 a cos 2 p 

n sin p 
2 a cos 2 p 


[q-E 


COS 2 (j7t/Z — 0) 


— (sin 2 (j7r/Z — 0) + tan 2 </?) 3 / 2 


[/<?<»]■ 


Obviously, fg(p) has the same sign with ^ and fg(p) is monotonic in¬ 
creasing. Since lim^i fg(p) = Ci > 0, we see that for given 6 € (0, ^), 
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(p £ [0, §), 

sup U(ip, 9 ) = max{[/( 0 , 6 ), U((p, 0 )}, 

¥>e[o,£] 

and the conclusion follows. 


□ 


Lemma 3.7. For all n £ N + , C n < -(logn + 7), then U( 0, f-) = < 


2an 


(log n + 7). 


Proof. Let 7 be the Euler-Mascheroni constant, that is 
7 = lim (1 + - 4-b — — log n). 

n—)-oo 2 n 

Let a n = 1 + 5 H-b 7 — log n and b n = 1 + \ J \ -b 731 — log n, then 

1 n + 1 

®n+l — 7 I 7 lo§ 7 ^ 0, 


n + 1 


n 


which implies 
Thus we get 


1 n + 1 

K+i -K = -log-> 0, 

n n 


b n < 7 < an- 


n— 1 


V - < logn + 7. 


3 = 1 

We notice that, for 0 < x < 1 , we have 

sin 7ix 

7T < 


x(l — x) 


< 4 , 


which implies 


then 


j 1 1 n 1 1 . 

CSC—7T <-:-= —-1-), 

71 ^h^-n) 77 J n ~ J 


^ n— 1 . 71—1 

^ 1 v j n v-^ 1 n 

c n = 7> cse-irC-) tC -(logn+ 7), 
2 n 7r 7 7T 

7=1 7=1‘ 

and the conclusion immediately follows from ( 13 . 31 ) . 


□ 


Remark 3.8. There is an asymptotic expansion of C n in Lemma 1 of m 
for large n: 


C n 


n 


7r 


2 n 


, _ All. ^ V—^ 

(7 + log—) + 2^ 
7r z —' 


( -l) fc (2 2fc -l - l)S 2 2 fc 7T : 


2 _ 2 fc —1 


fe>l 


{ 2 k){ 2 k)\ 


n 


1 

2fc—1 ’ 


where B2k stands for the Bernoulli numbers. 









16 


ZHIQIANG WANG AND SHIQING ZHANG 


Lemma 3.9. For n = 21 > 8 , we have U(t%, 0) < ^ (log ri + 7 ). 
Proof. From (I3.3p . we have 


l 


n j=l 


n 


1 -1 


—{2Q + \^[sin 2 (^) + tan 2 —] 2 + cot —} 
— ^' / T7. r?. 


4a cos : 


i=i 


-r 1 

n 


7 r. 
n J 


< 


n , 7 T, n , 4Ci 1 . 

4a cos - 1 COt n = 4a cos^ sin"^" 

n n n 


n 4 Ci 1 

< ~—(- 2 1 - 3 “)) 

“ — TT 13 7 

2n 2 n Qn 6 

Thus by Lemma 13.71 we have 


ml 0) , , 1 °gj+7 1 1 

n J < 2avr 1 _ jeL + 2 1 _ jeL 
2 rF 1 IFF 


n 

2air 


[(1 + 


2n 2 


1 - 


, n ,1 1 

^-Xlogj+7) + 2T— 

2n 2 6n 2 


n 2 n n TT n 11 

2^ [log 2 +7 + 2rfT^< log 2 + 7) + 2rr^- 1 ’ 


6n 2 


since n > 8 and 7 w 0.57721566490153286, 


U( l' 0) K h [ '° gn + 1 - ‘° g!2 + l28^ ( ‘° g4 + 7) + 


6x64 


< 


n 

2air 


(log n + 7 ). 


□ 


Then the estimate (|3.2I) holds from Lemma 13.6113.71 and 13.91 


Remark 3.10. To finish the proof of Remark 11.71 we only need to prove 
A(u) < B for n = 8 and s = h = 2. Throughout this remark, we keep in 
mind that n = 8 and B = 42 x 2 1 / 3 7 t 2 / 3 T 1 / 3 . It is direct computation that 


C/(|,0)<t/(0,|) = -C' 8 
o o a 


2 \ J 7T 
- > CSC — < 
a ^ 8 

3 =1 


24 

a 
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then by Lemma 13771 the potential of the the test loop U(u) < Similar 
to the proof of Proposition 13.41 we have 

/' T f T 36a 2 7r 2 24 

*4(«) = / < / [^775 - 1 - ]dt, 

Jo J o 1 ~ a 

we choose a = 3 _ 1 / 3 7 r~ 2 / 3 T 2 / 3 , then *4(fi) < 36 x 3 1 / 3 7 r 2 / 3 T 1 ' /3 < £>. 

3.2. Partial Collision. In this section we prove the following theorem and 
the idea is mainly from Fusco et al. [9]. 

Theorem 3.11. A minimizer u* of A\/c s is free of partial collisions. 

The proof is by contradiction. Suppose u* is a minimizer of A\k, s with 
partial collisions at time t c £ [0,T], then we prove that there is some local 
deformation which has lower action and so the minimizer is collision free. 
With the following lemma, we can only consider u* with partial collisions 
at time t c and without collisions in the time interval [t c — r, t c ) U(f c , t c + r]. 

Lemma 3.12. Suppose u * is a minimizer of A\jc s with partial collisions at 
time t c , then the collision is isolated. 

Proof. A collision is called isolated at t c means it is an isolated point in the 
set of collision times and the lemma is just the Corollary 5.12 in [8j, so we 
omit the proof. □ 

As mentioned in Remark 11.41 it is enough to consider t c £ I = [0, £]. 
Moreover we see that u* £ A s has the symmetry (11.21) . it is obvious that the 
collision must happen at F \ 0. So in the following, we only discuss partial 
collisions in two situations: 1 . colliding at ^ 3 -axis and it is two regular 
/-polygonal collisions (Figure 3); 2. colliding in the ^i^-plane and it is l 
binary collisions (Figure 4). 

We notice that u* £ IC S implies that there are both symmetry and 
topological constraints on uo(t) at time t = 0, For example, for e small, 
Vf £ (0, e), uo(t) = Riuo(-t) and uo(0) € P 0 - . This doesn’t allow general 
perturbations since uq(0) is not in the whole plane Pq, for this reason we can 
not use Marchal’s idea of averaging the action over sphere or its extension 
in [8j for averaging over suitable circles (rotating circle property). But for 
t G ( 0 , t^), such technique works, so we only put our focus on the partial 
collision at time t = 0 (and the case t = X is similar). 



18 


ZHIQIANG WANG AND SHIQING ZHANG 


£3 



£3 



6 


Let k C n be the colliding cluster, and q(t ) = (qj(t))j G k are the trajec¬ 
tories of the colliding cluster. Denote the partial Lagrangian £k(o , (^)) = 
it'k + Lie where A'k = ^2j£k^ m j\Qj\ 2 is the partial kinetic and U k = 
S?; jek i<j ^Xi~-x~\ ^ ie partial potential function . Thanks to the blow-up 

technique (Lemma 13.13llTT5jl . it is enough to consider a parabolic collision- 

2 

ejection orbit q(t ) = (Kf)ss instead of q(t), where s is a normalized central 
configurations. Also since the “blow-up” sends all other bodies not con¬ 
cerned by the collision cluster k to infinity, we can only do deformation 
inside k (for more detail, see Section 7 in [8]). 

Lemma 3.13. (Proposition 6.25 in [8]j Let k be a colliding cluster, I & = 
/Ljek m i ( il : Then there is k > 0 such that the following asymptotic estimate 
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hold: 

4 ~ (Kt)% 

■ 4 i 

4 ~ ) 3 

■ ■ 4 o ^2 

4 ~ rfi' («*) 3 • 
y 

and therefore 

1 •• 2 o -2 

/i k ~ £4 ~ ^4 ~ q k ( Kt ) 3 • 

_ _1 

Let s = 4 2 Q be the normalized configuration, we have 

Lemma 3.14. (Proposition 6.32 in 0 ; For every converging sequence s(tj ) 
of normalized configuration. The limit s = limj_>. 0 0 s(tj) is a central config¬ 
uration. 


We say that q is a (right) blow-up of the solution x(t) in 0, if 

q = (Kt)3s. 

For every A > 0 consider the path x x defined by 

x x (t) = X~3x(Xt) 
for every t € [0, A -1 e], we have 


Lemma 3.15. (Proposition 7.9 in 0 ; Let x(t) be a solution in (0, e), with 
an isolated collision in t = 0. Let k C n be a colliding cluster and q a (right) 
blow-up of x(t) with respect to k in 0. Let r £ (0, e) and let ip be a variation 
of the particles in k which is C 1 in a neighborhood of t, defined and centered 
for every t € [0, t\. Then there exists a sequence ip n € l4([0, r], M 3 ) fc of 
centered paths converging uniformly to 0 in [0, r] and with support in [0, t\, 
has the following property: 


lim 

n—> oo 


[ T [C(i 

Jo 


+ + if n ) - C(x Xn )\dt = / [C(q + ip) - C(q)\dt. 


Furthermore, our deformation is based on the solution of the two-body 
problem([2] [5]). Suppose u : M —» M 2 is the ejection-collision solution of 

f3u 


uj = — - 


M 


P > o, 


i.e. u(±t) = (|/3) 1 /4 2 / 3 n ± , t > 0, where n ± = lim^ 0 + ■ Let 

6d = arccos(n + • n~) and 0* = 2ir — 9d- If 0^ € (0,7r], given e > 0, there are 
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W) 


Figure 5 . 


exactly two Keplerian arcs u'■ [—e, e] —)• M 2 and Wj : [—e, e] —>• M 2 which 
connect uj(—t) to u(r) in the time interval [—r, r] and satisfy 


Wd((-e,e)) C {ain + a 2 n + , a, > 0}, 

Wj((—e, e)) C spanjn - , n + } \ {ain - + a 2 n + , a* > 0 }. 


The arcs and Ui are called the direct and indirect Keplerian arc. When 
6d = 0, the indirect arc does not exist and u>d((— e, e)) C {an + ,a > 0}. 

Lemma 3.16. The following inequalities hold: 

(1) A(u d ) < A(u |[_ e , e] ), Vn ± , 

(2) A(ui) < , Vn 1 * 1 such that n + ■ n~ < 1. 

where 



Proof. There is a detailed proof in Proposition 5.7 of [9J. 


□ 


3.2.1. colliding in ^-axis at time t = 0. In this situation, uo(0) € £3 fl Pq 
and since it is a partial collision, we must have uo(0) • e 3 < 0. We can do 
perturbation in the whole plane £3 = uq( 0) ■ e 3 . 
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The collision set at t = 0 is q = {qo, Qi, ■ ■ ■, Qi- 1} and it is a collision of 
regular /-polygon, so the partial functional 

[ C k (q)dt = l [ (^|<?o| 2 + Ar)dt, 

Jo Jo z I Qo I 

where a = \ csc(jn/l). Let q be the right blow-up of q, then % is the 
collision-ejection in the plane £3 = qo(0) • 63. By Lemma T3.161 there exists 
a Keplerian direct arc oj(t) satisfying J Q r ( 4 |u;o| 2 + ]^)dt < fo (II®1 2 + 
1^1 )dt and if we let q 0 (t) = (u(t),q 0 (t) • e 3 ), we have q 0 (r) = q 0 (r). By the 
symmetric condition (|1.3lh we must have %(0) € P 0 _ \^3. Consequently, the 
perturbation q = {%, Ri%, ■ ■ ■, Ri-iqo} satisfies / Q r C k (q)dt < / Q r C k (q)dt. 
By Lemma l3.13H3.15l there is no partial collisions in ^3-axis. 

3.2.2. colliding in the t^i^-plane at time t = 0. The partial collision is the 
union of l binary collisions. This situation is more complicated since we 
have the topological constraint uq( 0) ■ e^ < 0 which implies that co(t) should 
be the indirect arc when 0 < 0 < f (Figure 5). But the indirect arc does not 
exist when 9 = 0, so similar to the proof of Section [3.2.11 there is no partial 
collisions unless n ± = lim^ 0 ± = es ( see Fi S ure 4 )- 

Remark 3.17. Such collision described in the above is called the collision 
of type (=t) in Definition 5.1 of [9j. Let u* E K. be a minimize!' of the action 
A\k. and assume that u * has a partial collision at time t c . Let r be the axis 
on which the collision of the generating particle takes place and n + ,n~ be 
the unit vectors associated to the collision. Then we say that the collision 
is of type (=l) if 

(1) n + = n~, 

(2) The plane generated by r, n = is fixed by some reflection R € G. 
Lemma 3.18. (Corollary 5.1 in w 

Let 10 : (0, t) —» (—t, 0) —>• M 3 ) be a maximal ejection (collision) 

solution to the equation 

(o a\ ■■ (-^r ~~ 4 ) w T T ( \ 

(3 ' 6) " = ° |(JWM3 + Vl{w) - 

and let n = lim,_ (0+ be the unit vector or- 

thogonal to r. Assume that the plane 7iy jn , generated by r,n, is fixed by 
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some reflection R. Then 

£ 7 T rtn , Vi £ (0, t), (or Vi € (—i, 0)). 

In \S. (Xi R n denotes the rotation of it around the axis r, V\ is a smooth func¬ 
tion defined in an open set It C M 3 containing r \0 and a£l. Moreover V\ 
satisfies the symmetry condition Vi(Rlo) = RV\(uj), where R is a reflection 
such that Rr = r. 

By the above arguments, the collision of type is the only collision 
at time t = 0. Since there is no collision in (0, ^), we can apply this lemma 
to our context and let i = then uo(t) £ Pq, Vi £ (0, ^), which is a 
contradiction with uo(^) Ps i n our assumption, thus there is no partial 
collision at time t = 0. And the case for t = ^ is similar, so we have finished 
the proof. 

Remark 3.19. Lemma l3.18l is actually right in the view of physics. If there 
is some plane ir such that the particle force ^ £ 7r when u £ tt. Then 
w(0),w(0) £ 7r must imply that u(t) £ ir for all i > 0. 
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